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THE POPULATION OF M-TYPE SUPERGIANTS IN THE 



^ ; STARBURST CLUSTER STEPHENSON 2 

Ignacio Negueruela^, Carlos Gonzalez- Fernandez ^ , Ricardo Dorda^, 
' Amparo Marco^ and J. Simon Clark^ 

, Abstract. The open cluster Stephenson 2 contains the largest collection 

. ' of red supergiants known in the Galaxy, and at present is the second 

i-C ' most massive young cluster known in the Milky Way. We have obtained 

multi-epoch, intermediate-resolution spectra around the Can triplet 
^ ' for more than 30 red supergiants in Stephenson 2 and its surroundings. 

$-H , We find a clear separation between a majority of RSGs having spectral 

lyj ■ types M0-M2 and the brightest members in the NIR, which have very 

' late spectral types and show strong evidence for heavy mass loss. The 

distribution of spectral types is similar to that of RSGs in other clusters, 
such as NGC 7419, or associations, like Per OBI. The cluster data 
^ I strongly support the idea that heavy mass loss and maser emission 

■ is preferentially associated with late-M spectral types, suggesting that 

fT^ ' they represent an evolutionary phase. 

oo " 

ff^ 1 Introduction 

O ' 

' Red supergiants (RSGs) are very difficult to study in isolation. Accurate stellar 
. parameters depend on a good estimation of their distances, while comparison of 
^ ' these parameters to evolutionary models requires knowledge of their ages. Such 
\^ . information is not available for isolated RSGs. 

' Observing extragalactic populations (e.g., [Levesque fc Massey 2012] ) removes 
: the distance problem and provides large, statistically-significant samples, but may 
introduce other complications (fainter targets, spatially unresolved systems, con- 
tamination by bright AGB stars, possible metallicity gradients, mixed ages, ...). 

For Galactic stars, distances and ages are inferred from membership in open 
clusters or OB associations ( [Levesque et al. 2005[ ) . Unfortunately, given the short 
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duration of the RSG phase, typical open clusters in the solar neighbourhood, with 
total initial masses ~ 10'^ Mq, do not contain more than one or two RSGs. Indeed, 
simulations of open clusters using modern stellar evolution models and a standard 
IMF suggest that cluster masses ~ 10'' Mq are needed to have a good statistical 
chance of finding 2-4 RSGs (jClark et al. 2009a|) . Therefore population analysis 



must be done by combining large cluster samples (e.g., Eggenberger et al. 2002| 



Levesque et al. 2005) and thus mixing RSGs of different ages and perhaps even 



somewhat different chemical compositions. 

Mixed populations can be used to derive the current properties of the RSGs, but 
provide only weak constraints on their parent populations. Stronger constraints 
on evolutionary models could be obtained by analysing moderately-sized single- 
age populations of RSGs. The Perseus Arm contains some interesting examples of 



such populations (see Humphreys 19701. The Per OBI association includes almost 
20 RSGs, but its exact boundaries are not certain - not only its spatial extent, but 
also its actual membership (cf. IWalborn 2002^ . The core region, centred on the 
double cluster h and x Persei, contains ~ 10 RSGs for a total mass ~ 20 000AfQ, 
though again measuring the total mass over an area ~ 1° across is a difiicult 
process (jGurrie et al. 2010)) . 

The best example of a well-defined RSG population is provided by NGC 7419. 
This compact open cluster, with a mass > 5 x 10"^ Mq contains 5 RSGs at an age 
14 ±2 Myr ( [Marco fc Negueruela 2013 ). One of its most remarkable features is the 



presence of the very late supergiant MY Gep. This object, which presents strong 
evidence for heavy mass loss ( [Fawley 1974|Verheyen et al. 2012] ), seems to have 
been spectroscopically stable at spectral type M7.5I for the last ~ 50 yr. Such 
a late spectral type is very unusual for an RSG. Moreover, the few other RSGs 
observed with very late (>M4) spectral types tend to spend most of the time at 
earlier spectral types, and go to late-M types only during brief excursions (e.g., 
VX Sgr, which displays excursions up to MIO; see [Wing 2009[ ). 

In an analysis of the Milky Way RSG population, Levesque et al. (2005; see 
also Levesque, in these proceedings) found M-type supergiants with spectral types 
between MO and M5, with a clear peak around spectral type M2. This result is 
based on a sample of RSGs in OB associations. There is some hint that the RSGs in 
open clusters have a slightly different distribution of spectral types, preferentially 
having spectral types in the M0-M2 range, with a small fraction of them having 
later spectral types and indications of heavy mass loss. For example, in NGC 3766, 
the two RSG members have spectral type MO lb ( [Humphreys 1978] ). In x Per, 
FZ Per, B P -H56°595 and HP 14580 are MO-lIab, AP Per is M2.5Iab and RS Per 
is M4Iab ( [Humphreys T970| ). In NGG 7419, four RSGs are M0-2Iab and MY Cep 
is much later. Unfortunately, the number of open clusters with more than one 
RSG is small, and all these clusters contain relatively small populations of RSGs. 
Once again, coming to any conclusion implies combining data from many different 
clusters and ignoring potentially complicating factors, such as an age dependence. 
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2 Stephenson 2 

The advent of infrared detectors has allowed the exploration of the inner Milky 
Way, leading to the discovery of many sites of intense star formation and many 
massive young star clusters. In the last few years, several clusters containing many 
RSGs have been found towards the base of the Scutum Arm ( |Figer et al. 2006| 
IClark et al. 2009bl) . The open cluster Stephenson 2, though already known and 



identified as a candidate massive star cluster ( Stephenson 1990[ ), turned out to 
contain the highest number of RSGs (jPavies et al. 2007p . This cluster lies behind 
a dust layer (LDN 515) and is thus affected by heavy extinction {Ay « 11 mag). 

Davies et al. (2007) obtained A'-band spectra of sources within a circle of radius 
7', finding a large number of RSGs with similar radial velocities. By assuming that 
stars with radial velocities within ilOkms""'^ of the average value were cluster 
members, they concluded that 26 RSGs belonged to the cluster. The average 
radial velocity wlsr — +109.3 ± 0.7 kms~^ implies a kinematic distance « 6 kpc. 
If the cluster is at that distance, its age is ^ 17 Myr, and its total mass must 
be > 5 X lO'' Mq. From the depth of the CO bandhead at 2.3 /xm, Davies et al. 
(2007) estimated luminosity and spectral type for all their stars, as this feature is 
sensitive to both parameters. 

The reddening to Stephenson 2, though high, still allows observation of its 
RSGs in the region around the Can triplet. This spectral region provides much 
more information than the K band. Classification criteria for RSGs have been 



known since the photographic plate era (Keenan & Hynek 1945 Sharpless 1956) 



and have been extended to digital spectra ( [Carquillat et al. 1997[ ), though most 
of them are only valid for spectral types up to M3. Later types are not so well 
characterised, but criteria have been defined (e.g.. lSolf 19781 see also Dorda et al. in 
these proceedings). Observing a wide area around the cluster with the multi-object 
fibre spectrograph AutoFib2+WYFFOS (AF2) mounted on the Prime Focus of 
the 4.2 m William Herschel Telescope (WHT; La Palma, Spain), Negueruela et 
al. (2012) found many more RSGs with radial velocities similar to the cluster 
average. They concluded that the spatial extent of Stephenson 2 cannot be clearly 
defined. A compact core, with radius < 2.5' and containing about 20 RSGs, is 
surrounded by an extended association that merges into a general background with 
an important over-density of RSGs. 



3 Observations 

We have complemented the observations reported in Negueruela et al. (2012) with 
new observations at higher resolution. Spectra of most RSGs in Stephenson 2 were 
obtained during three runs in 2010, 2011 and 2012 with the fibre-fed dual-beam 
AAOmega spectrograph mounted on the 3.9 m Anglo-Australian Telescope (AAT) 
at the Australian Astronomical Observatory. The instrument was operated with 
the Two Degree Field ("2dF") multi-object system as front-end. We used the red 
arm, suited with a low-fringing CCD, equipped with the 1700D grating, blazed at 
lOOOOA, and specifically designed to observe the Can triplet region. This grating 
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provides a resolving power i? = 10 000 over slightly more than 400A. The exact 
wavelength range observed depends on the position of the target in the 2dF field. 

Further observations were obtained in June 2012 with the red arm of the long- 
slit ISIS double-beam spectrograph, also mounted on the WHT. We used grating 
R1200R centred on 8700A. The instrument was fitted with the Red+ CCD, which 
has very low fringing in this range and covered the 8400 -8950 A range with a 
nominal dispersion of 0.26A/pixel. We used a 0.9" slit, resulting in a resolving 
power 16 000. 



4 Red supergiants in Stephenson 2 

Of the 28 RSGs listed by Davies et al. (2007), we have observed 24. One supergiant 
(D49) is very heavily reddened, and is not detectable in the / band down to the 
DENIS limiting magnitude {i « 18 mag). The two faintest RSG members (D52 
and D72) are not detected in DENIS either (perhaps due to crowding). Finally, 
D23, which is an outlier in terms of spatial location and radial velocity, is too faint 
for a spectrum in this spectral region. 

Our data fully confirm the nature of all the other supergiants, validating the 
results obtained by Davies et al. (2007) with the CO bandhead technique. The 
spectral types derived from the X-band spectra, however, are only approximate. 
The CO-bandhead calibration results in a small average dispersion, but a given 



star may easily be wrong by 3 subtypes (Negueruela et al. 20121 



Two of the brightest supergiants, Dl and D4, have CO-bandhead radial veloc- 
ities not consistent with cluster membership, but not far from the cluster average 
(|Davies et al. 2007j) . Our radial velocities, however, indicate membership for Dl, 
the brightest star in the field (in the K band). The radial velocity of D4 is still 
outside the range adopted for cluster members, but there are some reasons to think 
that this star belongs to the cluster (see below). 

Assuming that the very high extinction to D49 is interstellar, Davies et al. 
(2007) calculated that this is the intrinsically brightest member of the cluster, 
slightly more luminous than D2. If Dl is indeed a member, it is much brighter 
than any of the two. With (J — Ks) ~ 7 mag, D49 is a very bright mid-IR source 
(unfortunately, it seems to be saturated in both the GLIMPSE and WISE cata- 
logues). It is thus very likely that this RSG is embedded in a thick envelope and 
actually self-absorbed. It displays SiO and II2O maser emission, generally consid- 



ered strong evidence of mass loss into an extended envelope (Deguchi et al. 2010 
IVerheyen et al. 2012D . 

Maser emission has also been detected from Dl and D2, but not from any 
other member ( [Deguchi et al. 20T0l [Verheyen et al. 2012D , suggesting that these 
objects have the highest mass loss rates amongst RSGs in Stephenson 2. Only 
the most luminous cluster members (in terms of bolometric luminosity) seem to 
present detectable maser emission. 

Our spectra (see Fig. 1) show that D2 has the latest spectral class in the 
cluster. It may be slightly variable, oscillating around M7I. Dl is very slightly 
earlier. D3, D4 and D5 have spectral type M5I, while D6 and D9 are M3.5Iab, 
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Fig. 1. Three RSGs in Steplienson 2 compared with the very late supergiant MY Cep, 
a rnernber of the open cluster NGC 7419. The shape of the spectra for stars in Stephen- 
son 2 is determined by the very strong reddening (Ay ~ 11 mag). The presence of VO 
bandheads at 8538, 8624A (shown on the spectrum of MY Cep) places D2 at a spectral 
type not earlier than M7. Given the smaller TiO bandhead at 886OA, Dl is slightly 
earlier (around M6). Other prominent TiO bandheads are marked on its spectrum. 
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showing a remarkable correlation between brightness and spectral type (note that 
the numbering system of lDavies et al. 20071 reflects brightness in Kg). The only 
exception is D8, which might not be a real member, as it does look like a very late- 
type giant. All stars fainter than Dll for which we have spectra (14 objects), have 
spectral types between MO I and M2.5 1, with most concentrating around Ml. 5 lab. 

The distribution of spectral types in Stephenson 2 is thus consistent with 
the information from other clusters, suggesting that the M I subtypes represent, 
to a first approximation, some sort of evolutionary sequence. Given the ob- 
served fractions, most stars should spend a substantial fraction of their time as 
RSGs at spectral type M1-M2 and then move to later spectral types, becoming 
brighter. Extreme mass loss, as reflected by maser emission or very bright mid- 
IR colours, seems to happen only at later spectral types. RSGs in the Per OBI 
association also show a similar correlation between spectral type and luminosity 
and between spectral type and mass loss, measured from the mid-infrared excess 
(jCohen fc Gaustad 1973)) . Field RSGs in the Carina arm display similar correla- 
tions ( [Humphreys &: Ney 1974[ ). 

This view would be consistent with recent results suggesting that all M-type 
RSGs have similar effective temperatures and that the strength of the TiO spec- 
trum, the main criterion used for spectral classification, is strongly correlated to 
luminosity rather than temperature (jPavies et al. 2013| and also in these pro- 
ceedings). The evolution towards later spectral types would then be caused by 
the heavier mass loss, which would result in the production of more extended 
molecular layers as the star's luminosity increases. 

The clusters of RSGs will be excellent laboratories to test these hypotheses, 
as they provide large populations of RSGs. In a given cluster, all RSGs are ex- 
pected to have the same age and thus similar initial masses (though the effects of 
high initial rotation may lead to some mixing of stars with different initial masses; 
see lEkstrom et al. 2012[) . The comparison between different clusters can help dis- 
entangle any age dependence of the average spectral types. For this, accurate 
distances and ages to the cluster still have to be determined. 
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